Gravitational collapse in Hofava-Lifshitz theory 
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We study gravitational collapse of a spherical fluid in Hoi^ava-Lifsliitz theory with the projectabil- 
ity condition and an arbitrary coupling constant A, where |A — 1| characterizes the deviation of the 
theory from general relativity in the infrared limit. The junction conditions across the surface of a 
collapsing star are derived under the (minimal) assumption that the junctions be mathematically 
meaningful in terms of generalized functions. When the collapsing star is made of a homogeneous 
and isotropic perfect fluid, and the external region is described by a stationary spacetime, the prob- 
lem reduces to the matching of six independent conditions. When the perfect fluid is pressureless (a 
dust fluid), it is found that such matching is possible only in the case A = 1. In this case, the external 
spacetime is described by the Schwarzschild (anti-) de Sitter solution written in Painleve-GuUstrand 
coordinates. Our treatment can be easily generalized to other versions of Hoi^ava-Lifshitz gravity 
or, more generally, to any model of a higher-order derivative gravity theory. 
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I. INTRODUCTION 

The study of gravitational collapse provides useful in- 
sights into the final fate of a massive star jl]. Within the 
framework of general relativity, the dynamical collapse of 
a homogeneous spherical dust cloud under its own grav- 
ity was first considered by Datt |7 and Oppenheimer and 
Snyder 0]. It was shown that it always leads to the for- 
mation of singularities. However, in a theory of quantum 
gravity, it is expected that the formation of singularities 
in a gravitational collapse is prevented by short-distance 
quantum effects. 

In this paper, we study this phenomenon (classically) 
in the context of the Hoi'ava theory of gravity [4^ . Since 
Hoi^ava's theory is motivated by the Lifshitz theory in 
solid state physics [S], it is often referred to as Hofava- 
Lifshitz (HL) theory. One of the essential ingredients 
of the theory is the inclusion of higher-dimensional spa- 
tial derivative operators which dominate in the ultravi- 
olet, making the theory power-counting renormalizable. 
The exclusion of higher-dimensional time derivative op- 
erators, on the other hand, guarantees that the theory 
is unitary (the problem of non-unitarity has plagued the 
quantization of gravity for a long time [B]). However, 
this asymmetrical treatment of the space and time vari- 
ables inevitably leads to the breaking of Lorentz symme- 
try. Although such a breaking is much less restricted by 
experiments in the gravitational sector than it is in the 
matter sector [7l [8] , the question of how to prevent the 
propagation of the Lorentz violations into the Standard 
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Model of particle physics remains challenging [S] . 

The breaking of Lorentz symmetry in the ultraviolet 
manifests itself in strongly anisotropic scalings of space 
and time. 



X -> £x, t^ £H. 



(1.1) 



In (3 -|- l)-dimensional spacetimes, HL theory is power- 
counting renormalizable provided that z > 3 ^ llOj . In 
this paper, we will assume that z — 5. At low energies, 
the theory is expected to flow to z = 1. In this limit the 
Lorentz invariance is "accidentally restored." 

The anisotropy between time and space men- 
tioned above is conveniently expressed in terms of 
the Arnowitt-Deser-Misner (ADM) decomposition [llj . 
N, iV', gij, {i, j = 1,2,3), which are, respectively, the 
lapse function, shift vector, and the three-dimensional 
metric defined on the leaves of constant time. The 
requirement that the foliation defined by these leaves 
be preserved by any gauge symmetry implies that the 
theory is covariant only under the action of the group 
Diff(Af , J^) of foliation-preserving diffeomorphisms. 



St ^ -fit), 6x'^-C{t,^). 



(1.2) 



As a consequence, an additional degree of freedom ap- 
pears in the gravitational sector - the spin-0 graviton. In 
order to be consistent with observations, this degree of 
freedom needs to decouple in the infrared (IR). Whether 
this decoupling takes place or not is still an open question 
|12| . Let us point out that the spin-0 mode is unstable 
in the Minkowski background in the original incarnation 
of HL theory [4]. If the projectability condition 



N = N{t) 



(1.3) 
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remains imposed, this instability persists in the gener- 
alization of HL theory in which additional higher-order 



operators are included [131 [Tl] (although in this case the 
de Sitter spacetime is stable [E]). Another potential 
complication of HL theory is that the theory becomes 
strongly coupled when energy is very low [16. However, 
as long as the theory is consistent with observations when 
the nonlinear effects are taken into account, this is not 
necessarily a problem (at least not classically) . A careful 
analysis shows that the theory is consistent with obser- 
vations in the vacuum spherically symmetry static case 
[TT] and in the cosmological setting [T8H20] . 

One way to overcome the above problems is to intro- 
duce an extra local U{1) symmetry, so that the total 
symmetry of the theory is enlarged to [H] 



C/(l) K Diff(M, T). 



(1.4) 



This is achieved by introducing a gauge field A and a 
Newtonian prepotential ip. One consequence of the U{1) 
symmetry is that the spin-0 gravitons are eliminated 
pni22j . As a result, all problems related to them, such 
as instability, strong coupling, and different propagation 
speeds in the gravitational sector, are resolved. The [/(I) 
symmetry was initially introduced in the case of A = 1, 
but the formalism was soon extended to the case of any A 
[23H2S]. In the presence of a U(l) symmetry, the consis- 
tency of HL theory with solar system tests and cosmol- 
ogy was systematically studied in pSH^ . In particular, 
it was shown in [29 that in order for the theory to be 
consistent with solar system tests, the gauge field A and 
the Newtonian prepotential (p must be part of the metric 
in the IR limit (this ensures that the line element ds^ is 
a scalar not only under Diff(M, J^) but also under the 
local U(l) symmetry). 

Another possibility is to give up the projectability con- 
dition ( |1.3[ ). This opens up for new operators to be in- 
cluded in the action, in particular, operators involving 
Gi = N,i/N [16]. In this way, all the problems men- 
tioned above can be avoided by properly choosing the 
coupling constants. However, since this leads to a theory 
with more than 70 independent coupling constants |30j . 
it makes the theory's predictive power questionable, al- 
though only five coupling constants are relevant in the 
infrared. 

A non-trivial generalization of the enlarged symme- 
try (1.4) to the nonprojectable case N = N{t,x) was 



recently presented in [22 132] • It was shown that, as in 
general relativity, the only degree of freedom of the model 
in the gravitational sector is the spin-2 massless graviton. 
Moreover, thanks to the elimination of the spin-0 gravi- 
tons, the physically viable range for the coupling con- 
stants is considerably enlarged, in comparison with the 
healthy extension [10 , where the extra U(l) symmetry 
is absent. Furthermore, the number of independent cou- 
pling constants is dramatically reduced from more than 
70 to 15. The consistency of the model with cosmol- 
ogy was recently established in [32ti34] . In the case with 
spherical symmetry, the model was shown to be consis- 
tent with solar system tests [35]. In contrast to the pro- 
jectable case, consistency can be achieved without tak- 



ing the gauge field A and Newtonian prepotential cp to 
be part of the metric. Finally, the duality between this 
version of HL theory and a non-relativistic quantum field 
theory was analyzed in |37j . and its embedding in string 
theory were constructed in [38] . 

In this paper, we study gravitational collapse of a 
spherical star with a finite radius in HL theory with the 
projectability condition and an arbitrary coupling con- 
stant A [2TH24] . In general relativity, there are two com- 
mon approaches for such studies. One approach relies 
on Israel's junction conditions [39 , which are essentially 
obtained by using the Gauss and Codazzi equations. An 
advantage of this method is that it applies well to the 
case where the coordinate systems inside and outside a 
collapsing body are different. Although Israel's method 
was initially developed only for non-null hypersurfaces, it 
was later generalized to the null hypersurface case [40] , 
For a recent review of this method, we refer to [41" and 
references therein. The other approach is originally due 
to Taub |42j and relies on distribution theory. In this ap- 
proach, although the coordinate systems inside and out- 
side the collapsing stars are taken to be the same, the 
null-hypersurface case can be easily included. Taub's ap- 
proach was widely used to study colliding gravitational 
waves and other related issues in general relativity ;43,. 

In this paper, we follow Taub's approach, as it turns 
out to be more convenient when dealing with higher-order 
derivatives. Moreover, in contrast to the case of general 
relativity, the foliation structure of the HL theory implies 
that the coordinate systems inside and outside of the 
collapsing star are unique. Thus, also from a technical 
point of view, Taub's method seems a natural choice for 
the study of a collapsing star of finite radius in HL theory. 

The paper is organized as follows: In Sec. II, we give 
a brief introduction to HL theory with the projectabil- 
ity condition and an arbitrary coupling constant A. In 
Sec. HI, we write down the field equations relevant for 
a spherical spacetime filled with a fluid. In Sec. IV, we 
generalize these equations to include the case where an 
infinitesimal thin matter shell appears on the surface of a 
collapsing star, and give explicitly all the necessary junc- 
tion conditions. This generalization is carried out under 
the only assumption that the junctions should be math- 
ematically meaningful in terms of generalized functions; 
therefore, in this sense the generalization is the most gen- 
eral. In Sec. V, we apply the junction conditions to the 
case where the collapsing star is made of a homogeneous 
and isotropic perfect fluid, while the external region is 
described by a stationary spacetime. When the perfect 
fluid is pressureless (a dust fluid), we find that matching 
is possible only in the case A = 1. In this case, the exter- 
nal spacetime is described by the Schwarzschild (anti-) 
de Sitter solution written in Painleve-Gullstrand coordi- 
nates |44j . In Sec. VI, we present our main results and 
conclusions. 

We would like to emphasize that our approach can be 
easily generalized to other versions of HL gravity or, more 
generally, to any model of a higher-order derivative grav- 



ity theory. 

II. GENERAL COVARIANT HL THEORY 

In this section, we give a brief introduction to HL the- 
ory with the projectabihty condition (1.3) and an arbi- 



trary coupHng constant A. For details, we refer readers to 
[24| . The fundamental variables are {N, iV*, gij, A, (p), 
which transform as 



;/£ 



SN = C'"^kN + Nf + Nf, 

5 gij = ViCi + VjCi + /3ij, 
5 A = C^^A + fA + fA, 
5^ = f^ + C^^ip, (2.1) 

under Diff(7\f, T), and as 

SaA 

5^N, 



NV^a, Sagij = 



-a, 



(2.2) 



under the local U(l) symmetry, where a is the generator 
of the U(l) symmetry. The total action is given by 



where g — det gij , and 

£k ^ K,jK'^~\K^, 



c^ 



(2.3) 



Ca 



A 



N 
Cx = (1-A) 



(2A,-i?), 



(vV)^ + 2-ft:vV 



(2.4) 



Here the coupling constant A^, which acts like a three- 
dimensional cosmological constant, has the dimension of 
(length)"^. The Ricci and Riemann terms all refer to the 
three-metric gij. Kij is the extrinsic curvature, and Qtj is 
the 3-dimensional "generalized" Einstein tensor defined 

by 



Ki 



2N 

iti'j 



{~g,j+S7,Nj+\7,N,), 



-gijR + Aggij. 



(2.5) 



Cm is the matter Lagrangian density and Cy denotes the 
potential part of the action given by 



Cv = C^9o+9iR+ -T^ 



{92R^ + g^R^JB}'] 



^ (.94^^ 



g^RR^jR^^+geR^RlR'^ 



-^ [g^RV^R + 38 (V.i?jfe) (V^i?^'^-)] , (2.6) 



which preserves the parity, where the coupling constants 
5s (s = 0, 1, 2, ... 8) are all dimensionless. The relativistic 
limit in the IR requires that 



51 = -1, C = 



1 



167rG' 
where G denotes the Newtonian constant. 



(2.7) 



Variation of the total action (2.3) with respect to the 



lapse function N{t) yields the Hamiltonian constraint 

(fx^ \Lk + Lv- vQ''V^^:jV - (1 - •^) (VV)^ 
^%ttG I d^x^j\ (2.8 



where 



j,^^{NCm) 



5N 



(2.9) 



Variation of the action with respect to the shift N''' 
yields the super-momentum constraint 



V^' 



^T^J - ^Q^J - (1 - A)5,j VV = 8^GJ„ (2.10) 



where the super-momentum tt*^ and matter current J* 
are defined as 



j^] ^ ^K'^ 



\Kg^\ r 



-N- 



.SC 



M 



5N, 



(2.11) 



Similarly, variations of the action with respect to Lp and 
A yield, respectively, 

G'' (k,, + V,Vj(^) + (1 - A) V2 (k + VV) 

= 87rGJ<^, (2.12) 

i?-2Ag = 87rGJA, (2.13) 



where 



, _ SCm . _J{NCm) 
J^= ^^, Ja = i- 



5A 



(2.14) 



On the other hand, variation with respect to gij leads to 
the dynamical equations 



^Y' 



a - ipg^^ - (1 - A).g'^ vV 



-2 (K^) 
- 2(1 - A) \{K + V^ip)V'V^ip + K'^V^ip 



+ (l-A)[2V('F^)-5^^Vfc^^ 



pV + p^J + p^^ 



SttGt'^, 



(2.15) 



where {K^Y' = K^'KJ, /(,,) = (/,, + /,,) /2, and 



^ A 






^ \aR'^~ (v'V^-g'^V^U 



E^, 






F^^ ^ ^ ^^ f^""^ ^JZasC'iFsf, (2.16) 

with ris = (2,0,-2,-2,-4,-4,-4,-4,-4). The 3- 
tensors {Fs)^^ and F^'^ „, are given by Eqs.(2.21)-(2.23) 
in [22], which, for the sake of the readers' convenience, 



are reproduced in Eqs.(A.l) and (A.2) of this paper. The 



stress 3-tensor r*^ is defined as 

,, _ 2 5 (V5^m) 



y/g hi. 



(2.17) 



The matter quantities (J*, J% J^, J a, t^-') satisfy 
the conservation laws 



d^x^ 



gkiT 



kl 



1 



Vff 



(Vff^*), 



-2ifJ, 



A 



{\^Ja 






0, (2.18) 



VV,, - ^ (V5^,),t - ^ (V/.iV, - V.iVfe) 

-^Vfcj'^- + J^V,<^-^V,A-0. (2.19) 

In general relativity, the four-dimensional energy- 
momentum tensor is defined as 



rpfXl> _ 



1 6 [V^^^Ai] 



Sgj^ 



(2.20) 



where ^,v — Q, 1, 2, 3, and 

g'i^^^-N^+N^N,, g'^^^N,. gf ^ g^r (2.21) 

Introducing the normal vector n^ to the hypersurface t = 
constant by 



one can decompose T^i/ as follows [45] : 



(2.22) 



(2.23) 



where hltJ is the projection operator defined by hliJ = 

g}iJ + n^n^. In the relativistic limit, one may make the 
following identification: 



U\ Ji, Tij) = {-2pH, Sj, S^j) . 



(2.24) 



III. SPHERICAL SPACETIMES FILLED WITH 
A FLUID 

Spherically symmetric static spacetimes in the frame- 
work of the HL theory with U(l) symmetry with or with- 
out the projectabilty condition are studied systematically 
in[26l[2a[29l[35l[36lll6l|47]. In particular, the ADM 
variables for spherically symmetric spacetimes with the 
projectability condition take the forms 



iV = l, AT* = ^;e^('''*)-'^('^'*\ 



(3.1) 



in the spherical coordinates x' = {r,9,(j)), where dft^ = 
dff^ + sin^ 9 d(jP . The diagonal case A^' = corresponds 
to fJ.{t,r) = — oo. On the other hand, using the C/(l) 
gauge freedom (2.2), without loss of generality, we set 



^ = 0, (3.2) 

which uniquely fixes the gauge. Then, we find that 






= Cx, F^^ 



0, 



2i^' .. .. 



Rij — ^ ..j^j 



-2i/ 



n 



c 



K 



Ca 



(1-A) 

^ [4. 

-t-A -I 
r 

2A\ 

3 



Tv' - (1 - e^" 
i>2 - 2vp!e^-'' 

i^-" - ^e^^^-") (2rp! + 1) 
|e-''^(l-2rj/')+Agr2-l 






Cv = J:/:^^) 



V ' 



(3.3) 



where a prime denotes the partial derivative with respect 
to r, ^ij = df6^ + siia^ OSfSf, and Cy s are given by 



Eq.(Al) in [48]. The Hamiltonian constraint (2.8) reads 
/ {Ck +Cv~ SttGJ*) e^r^dr = 0, (3.4) 



while the momentum constraint (2.10) yields 



(1 - A){e^-'' [r^ifi" + m" - m'^^') + 2{fi'r - 1)] 
-!>'r2|+2r(Ai/'e 



oM-i^ 



= -8TrGr^e-''+''v, 



(3.5) 



where 



It can also be shown that Eqs.(|2.12|) and (2.13) now read 



e- (A„r^-1) + 1 



,2./ ^^^2 



et'+''^' - e^^v 



2(u' -Agre^'']e''+'' 



+(1 - A){e2'' (-r2,>" + r2;>V' - 2r;>') 

+6''+"^ [r^ip'" + Sfi'^i" - fi'iy" - 3fi"iy' + fi'^ ~ 3i^' ^i 

+2/iV'2) + 2r(2^" - u" + 2[u' - /i')') + 2i/'] } 



2rv' 



e"" (Ky 



(3.6) 



1+1 



47rGr2e2''JA. (3.7) 



The dynamical equations (|2.15|), on the other hand, yield 

1 



(l-A)r 






1 



+e2M(^" + ^^'2_^V) 



2(/i' + Ai/') + (4A-3) 
Arf//' - mV) + (2A - 1)(2^' - v') 



g2p_2e''+^(A/i + z>) 
(3.8) 



+ -{i\ + l)r^i'^ 



e2M + -re^'^CA 



+ (A^+i(A + l);>2)re2'' 

- [(2A - l)/i + r^'(t' + A/i) + Ar(t'' + /i')] e"+'^ 

= - — (^ee + i^el + 87rGr%) , (3.9) 



where 



2 

-I 
r 



+ r(A' + v4i^') -ylfl-e^'^)]^^, (3.10) 

and Fij is given by Eq.(A4) in Appendix A. We define 
a fluid with pr = pe ds a. perfect fluid, which in general 
allows energy flow along a radial direction, i.e., v does 
not not necessarily vanish [49] 



The energy conservation law (|2.18| now reads 
dr e^ r 



Ph + {PH + 4pr) V 
+4 {v - vfi) - A^J^ - 2 ( j^ + vJj^ ] = 0, (3.11) 
while the momentum conservation (2.19[) yields 



2 1 

w/i' - (u' - pj.) {v - pr + Pe) - t^JaA' 



V + u(2;> — /i) 



= 0. 



(3.12) 



To relate the quantities J* , J* and r^j to the ones often 
used in general relativity, in addition to the normal vector 



n„ defined in Eq.(2.22|, we also introduce the spacelike 



unit vectors x^, 6^ and 0^ by 

0^ = r5% (ff.^rsmOS'l:. (3.13) 

In terms of these four unit vectors, the energy-momentum 
tensor for an anisotropic fluid can be written as 

Tf^iy = PHn^n^ + q{n^_,Xv + n^Xi-i) 

+PtXiiXu + Pe {df^di, + (j>f,(f>,y) , (3.14) 

where pn, q, Pr and pg denote, respectively, the energy 
density, heat flow along radial direction, radial, and tan- 
gential pressures, as measured by the observer with the 
four- velocity n^. This decomposition is consistent with 
the quantities J* and J* defined by 

pH = -\j\ v^e^q. (3.15) 

It should be noted that the definitions of the energy den- 
sity Ph, the radial pressure Pr and the heat flow q are 
different from the ones defined in a comoving frame in 
general relativity. We refer readers to Appendix B of 
liSl for details. 



IV. JUNCTION CONDITION ACROSS THE 
SURFACE OF A COLLAPSING SPHERE 

The surface E of a spherically symmetric collapsing 
star naturally divides the spacetime M into two regions, 
the internal and the external regions, denoted by M^ 
and M"*" respectively, as shown schematically in Fig. [l] 
The surface E = dM~ = — 9M+ is described by 



$(f,r) = 0. 



(4.1) 



where ^(t, r) = r^TZ{t). The spherical symmetry implies 
that the ADM variables on M take the form (3.1). 



1. Preliminaries 

We assume that the normal vector V$ to the hyper- 
surface E with components 

<^.,=6l-TZ6l (4.2) 

$'^ = e-2''(l - e^" - ne>'+'')5^ + (e^"" + 71)5^, 

is everywhere spacelike, i.e. 

$'^$ ;, = £-2^^ [1 - (e^ -f e'-Rf] > 0. (4.3) 

This is the case if 7t is small enough. We may then define 
the vector field N = V<i>/||V<i>||g in a neighborhood of E.^ 



^ It must not be confused with the lapse functi on, a s in the present 
case it is set to one, as one can see from Eq.l 3.1 1. 




N^ 



FIG. 1: The spacetime is divided into two regions, the internal 
M~ and external M+, where M~ = {a;'' : r < Tl{t)}, and 
M+ = {x'' : r > 7^(t)}. The surface r = 7^(i) is denoted by 



N has length one, i.e. NxN^ = 1, and the restriction of 
A^ to S is the outward pointing unit normal vector field 
on E. 

Let H{^) denote the Heaviside function defined by 



Hm 



1, $>0, 
h $ = 0, 
0, $ < 0, 



(4.4) 



and let (5($) denote the delta distribution with support 
on S. By definition, S{^) acts on a smooth test function 
ip G C°° {M) of compact support by 



(^($),(^)= / ^dE, 



(4.5) 



where dE = tjvVolg is the volume three-form induced 
by (7 on E and t^r denotes interior multiplication by N . 
The derivatives 5^"-'($), n > 1, of (5($) are defined in a 
standard way and the following relations are valid |50j : 



A^(")($) = ^ 



<5(<I>), 

5("+i)($), 



$(5(")($) == -n(5("-i)($). 



71 = 0,1,2,..., 
n = l,2,.... (4.6) 



If / is a function defined in a neighborhood of E, we 
define the distribution /i5("-^($) by letting it act on a test 
function ip by 



(/,S(")($),^) = (,S(")($),/^). 



(4.7) 



The product fd{^) is well defined whenever / is C" and 
it depends only on the restriction /js of / to E. More 
generally, the product /(5'^"^($) is well defined provided 
that / is C" and it depends only on the values of / and 
its partial derivatives of order < n evaluated on E. 
Let F be a distribution on M of the form 

n 

F = F+iJ($)+F-[l-iJ($)]+^F^™('=)(5('=)($), (4.8) 

/c=0 



where the F„'s are functions defined in a neighborhood 
of E while F+ and F~ are sufficiently smooth functions 
defined on M"*" and M~ respectively. We define the func- 
tion F^ on M by 

F" =F+H{<S>) + F-[l-H{'S?)], (4.9) 

and we define the jump [F]^ of F across E by 

[F]-{x) = F+{x)- F-{x), xeT,. (4.10) 

We will also need the fact that the equation F = is 
equivalent to the equations 



and 



F±(x)=0, x€M^, 






(4.11) 



-fe) 



fe=0 



(j-fc)! 9$J- 



0, 



< j < n, 



where ^^ acts on a function / by 

1 



5/ 



|V$| 



df-N, 



and, more generally, for any j > 1, 



(PI 



|V$| 



-Liyd f 



(4.12) 



(4.13) 



(4.14) 



A proof of this fact is given in Appendix B. 
For n = 3, the conditions in (4.12) are 



Qplm{3) 



plm(2) 



5$ 

g2plm{3) Qplm(2 

9$2 ^ a$ 



0, 



plr,i{l) 



= 0, 



/(fipl^ ^ d2plm(2) gplmjl) ^ ^^^^^^^ 



^. Distributional metric functions 



(4.15) 
= 0. 



The field equations (2.8) - (2.15) involve second-order 



derivatives of the metric coefficients with respect to t 
and sixth-order derivatives with respect to x*. Thus, one 
might require that the metric coefficients be C^ with re- 
spect to t and C^ with respect to a;*, where C" indicates 
that the first n derivatives exist and are continuous across 
the hypersurface $ = 0. However, this assumption elim- 
inates the important case of an infinitely thin shell of 
matter supported on E. Therefore, we will instead make 
the following weaker assumptions: (a) /j, and v are C^ in 
each of the regions Af + and AI^ up to the boundary E; 



(b) fji is C° across S; (c) v is C° with respect to t and 
C^ with respect to r across S. 

The above regularity assumptions ensure that the 
mathematically ill-defined products (5(<i>)^ and (5($)i/(<I>) 
do not appear in the field equations. Indeed, the terms 
in the field equations (3.4| - (3.12) that could lead to 
products of this type are 



A* 



MAi 



(4.16) 



Our assumptions imply that these terms may contain 
H{<^f but not (5($)2 or (5($)iJ($). 

In order to compute the derivatives of /i and i/, we note 
that 

V = v'^ = v+H{-^)^v-[l- H{<^)], (4.17) 

where the functions /i+ and v^ are C^ on M+, while 
the functions /^^ and v^ are C^ on M^ . Let Vs denote 
an open neighborhood of E. Let jx^ and i>"'" denote C^- 
extensions of //+ and i^"*" to M+ U V^- Let (1~ and i>~ 
denote C^-extensions of /i" and z^^ to M^ U l^s- Then 
the functions 



^i. 



A+ 



Ai 



(4.18) 



are defined on V^ and the following relations are valid on 
E whenever a + /3 < 5: 



A= M 



df^drP 

ga+p 



M 






(4.19) 



Since /i is C° across S, we find 






A*,tr = (Ai,tr)^+A,t'5(*), 
Ai.rr = (Ai.rr)^ +A,r'5(*). 



(4.20) 



Since j/ is C° across E, the derivatives of v^ and i^~ in 
any direction tangential to S must coincide when eval- 
uated on E. In particular, since the vector U defined 
by 

U^ = 5^+il5^, (4.21) 

is tangential to E (i.e. U^N\ — 0), we obtain 

= C/^Ka]" - [v,tV +n[v,rr- (4-22) 

But \v^r\~ = 0, because v is assumed to be C^ with 
respect to r. Thus [v.t]" — 0. It follows that v is in fact 



C^ across E. The same argument applied to v^t and v^r 
now implies that v is in fact C^ across E. We find 



/3) ^ 



V , 



\D 



\D 






D 



(4)N|D 



>(3)^($) 



1/(5))^ + 2j>('^),5($) + i}(3)<5'($), (4.23) 

where z^*^"^ = d^iy/dr^. We emphasize that the expres- 
sions on the right-hand sides of (4.20 1 and (4.23 ) are inde- 



pendent of the extensions used to define fi and ;> in (4.181, 



because the values of jl, v, and their partial derivatives of 
order < 5 are uniquely prescribed on E in view of (4.191. 



3. The junction conditions 

We will find the junction conditions across E by sub- 
stituting the expressions (4.20) and (4.23) for the deriva- 



tives of /i and v into the field equations (3.4) - (3.12) 



Suppose that the energy density pn = — 2J* has the 
form 

oo 

PH = (Pff)^ + E/^2"'"^'5<"^W: (4-24) 

n=0 

where it is understood that o nly fi nitely many of the 
Pj^ s are nonzero. Since, by (3.3), 



Ck = {CkY 



Cv = {CvY 



the Hamiltonian constraint (|3.4|) reads 
C 



r<n{t) 



{^K 



(£+ 



AiiGp^ 



^ e^r dr 



K 



lr>n{t) 

oo 

+47rG5](-l) 



£+ + 47rGp+) e^r^dr (4.25) 






r=n(t) 



I Im{n) y 2\ 



The left-hand sides of Eqs.(3.5), (3.6) and (3.7) have 



no supports on the hypersurface r = TZ{t). Thus, these 
equations remain unchanged in the regions M~^ and M~ , 
while on the hypersurface E they yield 



Jip — (Jif) , J A — {Ja) 



D 



(4.26) 



In fact, in order to avoid that the ill-defined product 
i?($)i5($) arises from the term J^A' in (3.12), we will 
assume that J a is C^ ■ 

The gauge field A has dimension [A\ — 4, so the action 
cannot contain terms like A" with n > 2, that is, it must 
be linear in A. We therefore assume that A has the form 

oo 

A{t,r) = A^ + ^A^"(")(5(")($). 



Variation w.r.t. 

lapse N{t) 

shift TV' 

gauge field A 
metric gij 



Name of equation 
Hamiltonian constraint 
Momentum constraint 



Dynamical equations 

Energy conservation law 

Momentum conservation law 



General version 



(2.101 



(2.121 



(2.131 



(2.81 



(2.151 



(2.181 



(2.191 



Spherically s ymm etric version 
(3.41 



(3.51 



(3.61 



(3.71 



(3.81 and (3.9 1 



(3.111 



(3.121 



TABLE I: A list of all field equations. 



Associated junction condition 



( 4.25 1 
( 4.26 1 
( 4.26 1 
(Q) 
4.271 and (4.28 
( 4.29 1 
(IsOl 



It follows that 

oo 

n=l 
A,rr = (A^rr)" + [2A^r + K"f^°^S{'S>) 

+[A + 2Af;"(") + A^,";(i)] (5'($) 

oo 



n=2 



Thus, 



n=l 



F,i = (^/.)^ + E^//'"'"^'5(")($), 



n=0 



where 



-2!^r„2 



[r^iA^rr)'' - '^,yiA,r)'' + HA^rY 



F, 



F, 



F, 



+ri^^rA^ - (1 - 62")^!^] , 
+r{A + .4^;^'(o)) + riy^rA^"'^°^ 

+ri/^^A^"(l) - (1 - e2'')yl^"(l)] , 

g-2z.r 2 fj^Im(n^2) _|_ 2^/m(n-l) 



ee 



A.Im{n) 
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+riy,^A^™(") - (1 - e2'^)A^"(")l , n > 2. 



From Eq.( A.4) we find that the functions {i^„}^^i con- 
tain no delta functions whereas 



{Fj)rr = {Fj) 



6e" 



+i>(3),5'($)]. 



Thus, (2.16) gives 



Frr = (i^..)^-(1657 + 6g8)^;>^'^5(<J>), 



-6i> 



Fee = {Feer~i8g7 + 3gs)- 



C 



Writing pr in the form 

oo 
n=0 



we find that Eq.(3.8) remains unchanged in the regions 



M+ and M , while on the hypersurface E it yields 



Y, \—[e-"'{l-2r,y')+Agr^-l]A 



Im{n) 






(4.27) 



Using (4.121, Eq.(4.27l can be rewritten as a hierarchy 



of scalar equations on S. 



Similarly, Eq.(3.9) remains unchanged in the regions 
M+ and M~, while on the hypersurface E it yields 

°° (-21/ 
n=0 '^ 



n=0 

2 
r 



oSi' 



r^/m(n) 
9G 



pAImin) ^ g^Q^2 Im(n) 



(5(")($)=0. 

(4.28) 



=0 

(i + Af™(o))(5($) 

OO 



n=0 



I- 



Note that 

PH,t = (PH,)'' + [p^"^-7^p«]^W 

OO 



and, by (4.26), 



w,t = {v,tr -nvS{<i>) 



(3.11) takes the form 



Thus, in view of (4.26), the energy conservation law 



Jdre'^r^{ {pH,r + [p^"/"^ - Ttp^] <5($) 



n=0 



+4(p,)^+4^p^"'("),5(")($) 



n=0 



'^.t 



+4((« O"" - 7^w<5($) - t;^A^,0 



that is. 



(/ +/ ]e^'r'^{pH,t + v,t{PH+'^Pr) 

\Jr<n{t) Jr>n{t)J 

-7tp^ + z.,,(p2"^°)+4p^™("))-47tt) 



r=n(t) 






"^^2/ -f"i(n) T^ /m(n-l) 

= 0. (4.29) 



The momentum conservation law (3.12) remains un- 



changed in Af + and AI while on the hypersurface E it 
yields 






(4.30) 



where we have used that 

OO 

p; = p,<5($) + Y, [(Pr"^"^),.--^^"^ ($) + p^™(")^("+i) ($)] . 

n=0 

This completes the general description of the junction 
conditions, which are summarized in Table 1. In what 
follows, we will consider some specific models of gravi- 
tational collapse for which the spacetime inside the col- 
lapsing sphere is described by the Friedman-Lemaitre- 
Robertson- Walker (FLRW) universe. 



V. GRAVITATIONAL COLLAPSE OF 

HOMOGENEOUS AND ISOTROPIC PERFECT 

FLUID 

In this section, we consider the gravitational collapse 
of a spherical cloud consisting of a homogeneous and 
isotropic perfect fluid, described by the FLRW universe, 



ds^ 



-dp +a^{t) 



dr"^ 
1 — fcr^ 



f^d^Q. 



where k — 0,±1. Letting r = a{t)r, t ^ t, the 
corresponding ADM variables take the form (3.1) with 
N" = 1, and 



1 / r2 

.~{t,r) = --ln(^l-fc-^ 



P {t,r) = In 



-d{t)r 



y/a^{t)-kr^J' 



(5.1) 



where d < for a collapsing cloud. For a perfect fluid, 
we assume that 



= Pr =P (t)^ w = 0. 



(5.2) 



We anticipate that the junction condition for v requires 
k = 0. Then, we flnd that 

zy"(t,r)=0, A*" (t, r) = In (-riJ), (fc = 0), (5.3) 

where H = d{t)/a{t), and that 



£-=£-=0, C^^2AgA~ 



(5.4) 
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It is easy to verify that the momentum constraint (3.5) is 



satisfied, whereas the equations (3.61 and (3.7 1 obtained 
by variation with respect to f and A respectively, reduce 
to 



SAgH + SnGJ- ^ 0, 
47rG J^ + Ag = 0. 



(5.5) 



It is convenient to consider the cases A = 1 and A 7^ 1 
separately. 

A. Gravitational Collapse with A = 1 

We first consider the case of A = 1. In this case, the 
static spherically symmetric exterior vacuum solution has 
the form [26] 



Since v = 0, we have F^- — —Kg^-, and the first dy- 



namical equation (3.8) reduces to the condition 

-a^A- + 2a^K„A- + 2(3A - l)ad + (3A - \)h^ 
r 

+ 2a^{8TTGp' -A) =0. 

If thi s co ndition is satisfied the second dynamical equa- 
tion (3.8 1 also holds provided that A'^ ~ rA^^^ = 0. On 



^+(r) 



1 



In 



2m-' 



1 

r 3 

2 



+ -Ar^ -2yl+(r) 



r 



A^(r')dr'\, 



tq 



v+ =0, 
for which we find that 



the omer hand, the momentum conservation law (3.12) 
reduces to J^A~ = 0. We conclude that the general so- 
lution when fc = is given by 



£+ =-A±- 2A, £+ = 2 A, C+ = 0, 



Jt = Jl 



Ph = 0^ 



(5.11) 



(5.12) 



J' 



3A,H 

•J A 



A„ 



SttG ^ 4t:G 

with A^ = A^ (t) being given by 

ij2 



(5.6) 



where m"*" , vq are constants and A'^ — A'^ (r) is a function 
of r only, yet to be determined. 

As mentioned previously, the condition that ly be C^ 
across S implies that fc = 0. We let the interior solution 



be of the form (5.3), and assume that the thin shell of 



AgA- + (3A - 1) 



-A = -8nGp. (5.7) 



matter separating the interior and exterior solutions is 
such that 



In the rest of this section, we consider only the case 



where Ag = 0. Then, Eq.(5.6) yields 



•J A ^ "^tfi ~ "' 



(5.8) 



A^A^ + A'"'^°^6{<i>), Ja 



PH=PH + p'H^'^mi 



0, 



(5.13) 



for which Eq.(l5Jl) shows that now A-{t) is an arbitrary where p^ - J^ - Pg - p 



0. 



function of i, and a{t) is given by 



(3A-1) - + 



A = -87rG'p. (5.9) 



Proposition V.l For the spacetime defined by [5^. 
(5.9), (5.11), the six junction conditions (4-25)- \4-3C 
reduce to the following six conditions: 



It is interesting to note that, since the Hamiltonian 
constraint is global, there is no analog of the Friedman 
equation in the current situation. This is in contrast 
to the case of HL cosmology [21], where a Friedman- 
like equation still exists, because of the homogeneity and 
isotropy of the whole universe.^ Although there is no 
analog of the Birkhoff theorem in HL theory, so that 
the spacetime outside the collapsing cloud can be either 
static or dynamical, we assume in this paper that the 
exterior solution is a static spherically symmetric vacuum 
spacetime. We also assume that the value of Ag is the 
same in the exterior and interior regions, i.e. 



-6i/^+2A + 47rGp^(i)) 



n(tf 



f-OO 

+4 / A\rdr + AnGp^^ 

J-R.{t) ' 

A{t,r) is continuous across S, 



/m(0)^2 



r=n(t) 



= 0, 



■R\'-^-H'^]{n- HTZ) - SnGp 



/m(0) 



(5.14) 
(5.15) 
(5.16) 

HU, 

(5.17) 



K 



A- 



0. 



(5.10) 



p^™(°)(i,7^(i)) = e-■'o^ 



^ Considering that homogeneity and isotropy are good approxi- 
mations for our observational universe, this global Hamiltonian 
constraint allows dust-like fluid to exist; this was first realized in 
|51l where it was considered as a candidate of dark matter. 



aJQ 



/o 



_ rt 2K(t) 
^ p Jo 

(It I 1 



dr 



p^"(")(o,7^(o)) 



His)nisrpH^,is) 
-n{s)p^{s)]ds 



(5.18) 



11 



rp + 2^^"^°^ = o« E. 



(5.19) Equation (5.261 then gives A|s = so that in fact A 



Moreover, the condition that ^ be continuous across S 
implies that 



is continuous across S, which proves (5.16). Equation 



^^ = 



A -3^2 



-nn - HH.fTz^ 



(5.20) 



(|5.22|) can now be written as 

A,r(e2'' ^ -^g^^ j^ 2i,^ + STrGpe"'"'] ^($) 

+ i(5'($) = 0. (5.27) 



Proof. For the spacetime defined by (5.12|-(5.13l, condi- 
tion (4.251 reduces to 



In view of (4.15) this yields 



f-K(t) 



{-&H{tf + 2K + ATiGp^[t)) / r'^dr 

Jo 

f-OO 

+4 / A\rdr + AttCp^ 
Jnu) 



pA^ + '^e'') + 2^,r + SttGp, 



lm(0) 



dA 



on S. (5.28) 



j™(0)^2 



r=7?,(t) 



0, 



Now observe that if a function f{t, r) is C" across E, 
then 



5/ 9/ 



on E. 



(5.29) 



which yields (5.14). Moreover, condition (|4.26| reduces 
immediately to (|5.15). 



Conditions (4.27) and (4.28) reduce to 



9$ dr 

Indeed, the continuity of / implies that the derivative 
of / in any direction tangential to E must vanish when 
evaluated on E; thus f t+Ti-f.r = on E. A computation 



F;^^™(")^($) + i^;^^"(i)<5'($) = 0, (5.21) using (Ipl), (ig, and Km now gives (Ip9) 



and 



On the other hand, since 



r(A,re^'' + 7t/i,re^)(5($) 



p+ = -iAr-2A+)e-'^^ 



-^Alm{0) 



F^;""^'"5{^) + -F^;"''-'>6'{^) 



AIm{l)^/, 



we find 



-F 



AIm{2)f-„ 



(5"($) + 8^Grpg™^°^(S($) = 0, (5.22) 



t^,r = ^{^r-2A+)e- 



(5.30) 



respectively, where we have used that 

Frr — (Frr) , Peg — {Fgg) . 



(5.23) 



Now 



Inserting the equations ( |5.29[ ) and (5.30) into ( |5.28[ ), 
we find 

Q(Ar - 2A+)e-2M _ l\ (^^^t^ + n^f^) + i^^ 
+8^Gp'™^"' = on E. 



F^ = 



+ A^™(°)y($) 

F^g = (F.1)^+E^«l'""^"^^^"'W' 



Since A^^ = A+ = A ^Tt ^, simplification yields (5.17) 
Condition (4.29) reduces to 



(5.24) 
(5.25) 



n(t) 

d_ 
dr 



e^r pfj ^dr + r 



PH,t - T^PH 



where 



=K(t) 



r^Ttp^"^"' 



r=-K(t) 

(5.31) 



That is. 



(Fgi)^ = r2(A..)^+r(4.)^, 
pA.imiD ^ r\A + 2A';^^'>y) + rA'"<°\ 



^{t)PHAt) Z-^^*^ 3 



r'^dr 



F, 



A,Im{2) 



^2^/m(0)_ 



Thus, equation ( 5.21[ ) can be written as 

[A + A^'"(°)] (5($) + ^^"(o),5'($) = 0. (5.26) 



Thus, by ( |4.15[ ), A^™(")|s = 0. Hence, yl^'"(o)(5($) = 
which gives 

= (A^™(o)j($))^^ = A^;"W,5($) + A^"(°)(5'($). 



a(0 JO 

+n{tf[p'^}'\t,n{t)) + n{t)p-^{t) 

+2n{t)n{t)p'j^'^°\t,n{t)) 

+n(tfTZ{t)PHT\t^nt))^0- 

Consequently, 

iI(i)p^,(i)7^(t)2 d r ,„^(o^ 



(5.32) 



rft 



PH^'it^nt)) 



-2'^^p'^^°\t,n{t)) + n{t)p^{t)^o. (5.33) 
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Solving this differential equation for p^ , we find 
(pSl). 



Condition (4.30) reduces to 



p+2p/™(o)^^(<j,) = 0. 



This yields (5.191 



Finally, the condition that fi be continuous across S 
can be written as 



nit) 



/ A+{r')dr' = H^'R^. 



(5.34) 



Since A is continuous acro ss S, we have A'^{7l{t)) = 
A~{t). Hence, multiplying (5.34) by TZ and then differ- 
entiating with respect to i, we find 



A7^^7^ - 2A77^ = 2HH_tTV + iH^Tl^Tl. 



Solving this equation for A j, we find (5.20). 



D 



The conditions (5.171 and (5.20) imply that 



A-3iJ2 



-n - HH^tTl 



nil + {2H'^ + 2H,t - A)7^ - 167rGp, 



7m(0) 



0. 



Solving this equation for TZ{t) we find the following equa- 
tion which expresses TZit) in terms of H{t) and the pres- 



sure Pg™' on the shell: 



TZ{t) = e-^ons)ds\j^^Q^ 



-WttG I e 

'0 



/m(0)/ 



Jo Hr)drPj^_^ns)l^^y ^^3^^ 



where I{t) is defined by 



I = 2H 



2gt A 
H H' 



(5.36) 



B. Dust Collapse with A = 1 

Suppose now that the perfect fluid in the interior re- 
gion consists of dust, i.e. 



Pr ^Pg = 0. 



Then, the condition (5.19) implies that 



lrn(0) 



0. 



(5.37) 



(5.38) 



Solving equation (5.9) for a(t) we find 
a(i) 



flQ cosh 3 



/3A/ 



ao{to~t)^/^ 



(t-to)), A^O, 
A = 0, 



(5.39) 



where ag and to are constants. In the following, let us 
consider the cases A 7^ and A = 0, separately. 



1. A>0 

In this case, substituting the expression for a(t) into 
(5.36) we obtain 



lit) 



tanh 



(io - 1) 



and then (5.35) yields 



n{t) =7^ocosh^ 



/3A 



(^0 - t) 



(5.40) 



where TZq is a constant. Condition (5.20) now implies 
that ^7 — 0, i.e. A~{t) = ^0 for some constant Aq. 



Then, by ( |5.16| ), A+{n{t)) = Aq. That is, A+{r) = Aq 
for all r such that r = 7?.(t) for some t. Hence, the form 
of (5.40 ) implies that v4+ = Aq for all (i, r) in the exterior 
region. This gives 



A(t, r) - Ao 



(5.41) 



Condition (5.14) now implies 

(-6iJ2 + 2A + A7:Gp^{t))n(t) 



PH^'\t.nit)) 



127rG 
A + 7rG[l + cosh(\/3A(io - t))]PH{t) 



67rGcosh3(^(<o-i)) 



(5.42) 



Subst itutin g this into condition (5.18), or its equivalent 



form (5.33), we infer that Pjj{t) satisfies 



12 



cosh =5 



3A, 



(io -^) i 4 cosh 3 



/3A 



[h ~ t) 



-7^o\/3Asinhf ^(to - t) ) J.p^^(i) = 0, 



PhW^Ph 



(0) 



(5.43) 



JO) 



where p\f is a constant. All the conditions of Proposition 
|V.1| are now satisfied. It only remains to consider the 
condition that p be continuous across S. This condition 
reduces to 







— + hn' 2Ao + 1^„(7^ - to) H^n' 



n 



6Aoro + 7^gA 



37^oCOsh^/3(y|A(^^_^)) 
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That is, the parameter tq is fixed by 

6771+ 



ro 



7^3A 



6Ao 



This imphes that 



+ 1, f^^^ 



A^\ 
3r ' 



(5.44) 



(5.45) 



Since all the field equations and junction conditions are 
now satisfied we have proved the following result. 

Proposition V.2 Hofava-Lifshitz gravity admits the 
following explicit solution when A = 1 and A > O.' 



+ 1, Mr2 


ir 


Y fi- =ln{-H{t)r 


^ = 0, Hit)^ ^tanh(^(to o), 


7^(i)=7^oCoshi('^(to-t)Y (^ 


Pr=pe=0, PhW^Ph^ A{t,r)^Ao, 


lm(0) . . , 

Ph ''■^ gjuen by ( 


5.42 


), 



(5.46) 



where to, TZq, Aq, and p\j are constants. 

For t < to the dust cloud is contracting. Ast /^ to, the 
radius of the dust sphere approaches its minimal value of 
TZ — TZo at t = ioi ^-nd the function e^ approaches zero: 



lim n{t) =7^o, 



lim e^^(*) 

tZ-to 



0, 



as shown schematically in Fig. [2j After the star col- 
lapses to this point, it is not clear how spacetime evo- 
lutes, be cause /x+ becomes unbounded as one can see 
from Eq.(5.46l, for which the extrinsic scalar K'^ , 



K+{r) = 



pUr) 



(5.47) 



also becomes unbounded, which indicates the existence 
of a scalar singularity at this point [55]. However, such 
a singularity is weak. In particular, the corresponding 
four-dimensional Ricci scalar remains finite, ^^'R = 4A. 
Thus, it is not clear whether the spacetime across this 
point is extendabl e or n ot. 



not both be positive 



/m(0) 



In addition, Eq.(5.42| shows that p^ and p^ can- 



-A7?.q/6 we can write /i+ in the form. 



b understand this, letting M — 



P 



1, f2M 

7: In 

2 V r 



Ar2 



(5.48) 



However, this is nothing but the Schwarzschild-de Sitter 
solution with mass M and a cosmological constant A, 
where M is negative. 




FIG. 2: The evolution of the surface of the collapsing star 
for A = 1 and A > 0, given by Eq.( |5.46^ . At the moment 
t = to, the star collapses to its minimal radius TZ{to) — TZo, at 
which the extrinsic curvature K^ becomes unbounded, while 
the four-dimensional Ricci scalar remains finite. 



2. A<0 

In this case, substituting the expression for a{t) into 
(5.36) we obtain 



m 



tan 



\t~to) 



and then (5.35) yields 

7^(t) = 7^o COS5 



-{t-to) 



(5.49) 



where 7?.o is another constant. Condition (5.20) now im- 
plies that A~^ = 0, i.e. A~{t) = Ao for some constant Aq. 
Then, by \b.l&\, A+{TZ{t)) = Aq. That is, A+{r) = Ao 



for all r such that r = TZ{t) for some t. We will assume 
that A+ = Ao for all (i, r) in the exterior region, i.e. 
A{t,r) — Ao- Condition (5.14) now implies 



p^"(°)(t,7e(t)) = 

|A|-7rG[l 



(-6iJ2 + 2A + 47^G/9^(t))7^(i) 



= 7^n 



127rG 

^i^/mit-to))]PHit) 

/3|A|, 



(5.50) 



6nGcos3{^^{t-to)) 



Subst itutin g this into condition (5.18), or its equivalent 



form (5.33), we infer that p^(i) satisfies 



PhW = Ph\ 



(5.51) 



JO) 



where p\j is a constant. All the conditions of Proposi- 
tion |VT] are now satisfied, while the condition that jjl be 
continuous across E reduces to 







2m^ 



1 



A7^2 - 2Ao + 7^Ao{n - ro) - H^TZ^ 



6m+ - 6yloro + 7^j]A 

37^oCOs2/3(Vp(i„i^^)) 
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Thus, the parameter tq is fixed by 
6to+ + 7^gA 



ro 



6Ao 



This imphes that 



+ 1 f2M |A| 2 



(5.52) 



(5.53) 



where M = |A|7?,q/6. Clearly, this corresponds to the 
Schwarzschild-anti-de Sitter solution. For /i+ to be real, 
we must assume that r < TZq. Similar to the last case, the 
extrinsic curvature K~^ at r = TZq becomes unbounded, 
while the four-dimensional Ricci scalar (■*)/? remains con- 
stant. Thus, in this case it is also not clear whether or 
not the spacetime is extendable cross r = TZq. 

In any case, all the field equations and junction condi- 
tions are now satisfied for r < TZq, and we have proved 
the following result. 

Proposition V.3 Hofava-Lifshitz gravity admits the 
following explicit solution when A = 1 and A < 0.' 



M+ = 


2 \3r 


(7^i]-r3)), 


P = ln(- 


-H{t)r) 


v = i 


), H{t) - 


= ^ 3 tan 


\ 2 ^* 


\ 

-to) , 
J 


n{t) 


= TZf) cos' 




m), 


(5 


Pr = 


P9 = 0, 


PHit)-pP, 


A{t,r)^ 


^0, 



^ is given by {5.50), 



Ph 



(0) 



where t^, TZq, Aq, and p^ are constants. 

The evolution of the surface of the collapsing star is 
illustrated in Fig. |3] The collapse starts at an initial 
time ti < to, and at time t = ts, the star collapses to 
a central singularity at which we have TZ{ts) = 0, where 
tg = io + 7r/^3|A|. Equation (5.501 shows that now 
both PiP and p^ can be positive, provided that |A| > 
27rGpP. 



3. A = 



In this case, substituting the expression (5.391 for a{t) 
into (5.36) we obtain 



lit) 



3(io - 1) ' 



and then (5.35) yields 

n{t) = TZoito - t)i , 



(5.55) 



where TZq is a constant. Condition (5.201 now implies 
that Al = 0, i.e. A~{t) = Aq for some constant Aq. 



r 




R(t) 




Rq 




v^ 










tv, 


ts 


t 



FIG. 3: The evolution of the surface of the collapsing star 
for A = 1 and A < 0, given by Eq.(5.54|. The star starts to 
collapse at a time t — ti > to. At the later time t — ts, at 
which TZits) = 0, the star collapses and a central singularity 
is formed. 



Then, by (5.16), A+{n{t)) = Aq. That is, A+{r) ^ Aq 



for all r such that r = TZ{t) for some t. Hence (5.55) 
implies that A^ = Aq for all (t, r) in the exterior region. 
Thus, in the present case we also have A{t,r) = Aq. 
Condition (5.14) now implies 



7m(0) 



pr'"it,nt)) = - 



{-6H^ + 4TTGp]j{t))TZ{t) 



12nG 



n. 



2-3G^(to-t)Vg(0 
' 9G7r(to - t)4/3 



Substituting this into condition (5.18), or its equivalent 
form (5.33), we infer that 



PH^t) 



(0) 

Ph^ 



(5.56) 



.(0) 



where p\j is a constant. All the conditions of Proposi- 
tion V.l are now satisfied, and the condition that p be 
continuous across E becomes 



= 



2m+ 

n 

„9m^ 



2Ao 



n 

9Aoro - 27^g 



Ao{n-ro)-H^n^ 



97^(i) 
Hence, the parameter tq is fixed to 



9m+ 



''o 



27^3 



9^0 



which implies that 
^ 1 



P 



■In 



0, N+^l, 



(5.57) 



(5.58) 



where r^ = 47?.o/9. This is nothing but is 

the Schwarzschild solution written in the Painleve- 
GuUstrand coordinates |44j. All the field equations and 
junction conditions are satisfied, so we have proved the 
following result. 
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K*= QO 



FIG. 4: The evolution of the surface of the collapsing star 
for A = 1 and A = 0, given by Eq. (5.591. At the moment 
t — ti < to, the star starts to collapse until the moment t = to, 
at which we have 7?. (to) — 0, whereby a central singularity is 
formed. 



Proposition V.4 Hofava-Lifshitz 
following explicit solution when A = 



M 



2 \ r 



M 



In I 



gravity admits the 
1 and A = 0.- 

-H{t)r), 



iy = 0, H{t) = --— 

6[tQ 

7^(t) = 7^o(to-t)^ 

Pr=Pe=0, Pffit) = 



ty 



(0) 



(5.59) 



A(t,r) = Ao, 



/m(o) ^ 2 - 3G'7r(to 
Ph =7^o- 



t?pf 



9G7r(to - t)^/3 



where tg, TZq, Aq, 



(0) 



Ph 



are constants. 



The evolution of the surface of the collapsing star is 
shown in Fig. |4] The star begins to collapse at the 
moment ti with a radius TZi[= TZ{ti)] until the moment 
t = io; at which we have TZ{to) = and a central sin- 
gularity is formed. The spacetime outside of the star is 
given by the Schwarzschild solution. Thus, as in GR, the 
Schwarzschild spacetime can be formed by the collapse 
of a homogenous and isotropic dust perfect fluid [T] . We 
note that p^ > for 



^0 " \ ^r< (0) <^<^o- 



C. Gravitational Collapse with A 7^ 1 

We now consider the case of A 7^ 1. For an exterior 
static spherically symme tric vacuum spacetime with A 7^ 
1 and Kg = 0, equation (3.7) implies that 



1, / 2B\ 

= n 1 , 

2 I r y 



(5.60) 



where J5 is a constant. On the other hand, for the interior 
FLRW region, we have 

v = - - ui I 1 - fe ^ , , , I . (5.61) 



1 / r 
= ^ t:: In 1 — fc „, , 

2 V a (*) 



Hence, the condition v 



+ - ,,- 



on S implies that = 



kTZ{t)^. Consequently, in order for a solution with Tl{t) ^ 
to exist, we must have fc = 0. The conditions that ly and 
ly^r be continuous across E then reduce to 2B/TZ{t) = 0. 
Thus, in order for a nontrivial solution to exist we must 
have k = B = 0. Thus, we have 



= 1^+ = 0, 



M 



In (-rH) 



(5.62) 



On the other hand, since A 7^ 1, the momentum con- 
straint (3.5) yields 



/i = /i (r) = In Cir 



(5.63) 



where Ci and C2 are constants. The field equations (3.6) 
- (3.9) are then satisfied provided that 



. + . ^ .+ ^'^1 3(1-3A)C2 + 2A 2 ,,,,, 
A+{r) = A+-^ + ^ p r^ (5.64) 

where Aq is a constant. It is interesting to note that 
this class of solutions was first found in 29J in the IR 
limit. However, since the restriction of the spacetime to 



the leaves t = constant is flat, we have Ri 



0, and 



the higher-order derivative terms of Rij vanish identi- 
cally, so it is also solutions of the full theory. Note also 
that the asymptotical-flatness condition requires Ci = 0. 
However, in the following we leave this possibility open. 
Then, we flnd that 



C+=3Cf{l-3X) 



GCl 



C+ = 2A, £+ = 0, 



0, 



Jl 



0, J+ =0, p+= 0. 



(5.65) 



In order for the integral over the exterior region in the 
Hamiltonian constraint (4.25) to converge, we also need 
to assume that 

3Ci2(l-3A) + 2A = 0. (5.66) 

Similar to the case with A = 1, the interior solution is 



still of the form (5.8), i.e. 



J^ = J-=0, A-=A-{t), 



with a being given by Eq.(5.9). In view of (5.4), we have 

.^-3(l-3A)i72, 
y- = 0. (5.67) 



C 



V 



0, 



c 



X 



0, 



C-, 



£y = 2A, £7 - 0, 



We assume that the thin shell of matter separating the 
interior and exterior solutions is such that 



Pr 



Pr 



/m(0) 



J^ = 4™(°)<5($), 



pr'^'sm, pH = PH+pr"'6m 



lm{0) , 



A = A^+A^'"(o)j($), Ja = 0, Ca^C^, (5.68) 
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again with p'^ = p~ — P~ (t) and pjj = pHit). For the 
sake of simphcity, in the following let us consider the 
above three cases separately. 



Proposition V.5 For the spacetime defined by (5.62)- 
(5.68), the six junction conditions ^^25)-(4.30) reduce 
to the following six conditions: 



(3(1 - 3X)H{tf + 2A + 4TTGp„{t)) 

^2 



-^^-4^Gp^"("V2 



nt)' 



r=n{t) 



3 

= 0, 



J-Hu; = 0, 



/m(0) ^ 

A{t^r) is continuous across S, 
3C|(1-2A) 3AC27^ SACiCa 



(5.69) 

(5.70) 
(5.71) 



2r5 



(■t 2K(t) , 



/m(0) 



(5.72) 



p^'"(°\t,7^(i)) ^e-^'^^S^'^^ p^"(°)(0,7^(0)) 



^(s)Pff(s) (is 



rp 



2p, 



/m(0) 



on E. 



(5.73) 
(5.74) 



Moreover, the condition that p he continuous across E 
can he written as 



n{t) 



Co 



Ci+H 



(5.75) 



Proof. For the spacetime defined by (5.62|-(5.68l, condi- 
tion (4.251 reduces to 



3(1 - 3A)d(i)2 



2A + 47rGp^(i) 



K(t) 



a{tr 

+ r Ucl{l - 3A) + ^ + 2K\ r'dr 
Jn(t) \ ^ / 

0, 



r dr 



, . ^ lm(0) 2 

-|-47rGp^ ^ 'r 



r=n{t) 



which, in view of (5.66), yields (5.69). Moreover, condi 



tion (4.261 reduces immediately to (5.70) 

The functions F and F^ are given by (|5.23|) - \b.2b\ 



also for A ^ 1. Hence, condition ( |4.27 ) implies that A is 
continuous across E just like in the case of A = 1. This 
proves (5.71[). Condition (4.28) then reduces to 



\r{p^re^^' +'Rp.^re^)6{<^) 

1 „AIm(Vl\ „. , ^ 1 

r 



-F, 
r 

'-F. 
r 



■4^™(0)5(ci,) + >,f "(^)<5'($) 



>5"{^) + 8^Grp;™^^^J($) = 0. 



That is, 

AA,r(e2^ + ile^') + 2A,r + 87rGpe"^°^j <5($) 

+ iy($) = 0. (5.76) 

In view of ( |4.15[ ) and ( |5.29[ ) this yields 

AA,r(e^^ + T^e'^) + i,,- + 87rGp^"^°^ =0 on E. (5.77) 
Since 



Ai.r = 



Gir3 - 2G2 

Gir4 + G2r' 



M,r = 



we find 



-3G2 



M,r = 



Gir4 + G2r' 



(5.78) 



Equations (5.77) and (5.78) imply 



Cir^ + C2r r) 
+87rGp^"*°^ = on E, 



which upon simplification yields (5.72 1 



Conditions (4.29) and (4.301 reduce to (5.73) and 



(5.74) as in the case of A = 1. 



Finally, the condition that p be continuous across E 
can be written as 



C{R.(t) 



C2 



nty 



-Hn{t), 



which yields (5.751 



D. Dust Collapse with A 7^ 1 

Suppose now that the perfect fiuid in the interi or re- 
gion consists of dust, i.e. p~ = 0. Condition (5.74) 
implies that 



/m(0) 



0. 



Solving equation (5.9) for a(t) we find 



aocosh^^ {uj{t-to)), At^O, A^^i, 
a{t)^{aoito-t)y^ A = 0, A^|, (5.79) 



a(t) arbitrary. 



A = 0, A 



where to 



3A 



2(3A-1) 



pr , and flg and ig ^^^ constants. 
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1. A/0, A/i 



Assuming that C2 7^ 0, condition (5.72) becomes 



Setting 7?.(t) = e^'-^i*u(t), the above equation reads, 

_ C2(1-2A) 3c,t 

2Au2 



which has the general solution, 



1/3 



where Cq^Iq are constants. Thus, we have 



7^(i) 



^^^ii^^-Coe-3C^^*V^', Ci^O, 



2ACi 



(^^^^^^(t-M 



1/3 



Ci = 0. 



On the other hand, (5.75) implies that 

C2 



nt) = [-^ . 2 



Ci + |a;tanh(w(t- ip)) 



(5.81) 
(5.82) 



The expressions (5.81) and (5.82) for TZ{t) are not com- 



patible; thus there are no dust solutions in this case. 



2. A = 0, A / i 



Condition (5.72 1 reduces to (|5.80[) with solution given 



by (5.81). On the other hand, (15. 751) implies that 




7^(t) 



Co 



Ci 



3(to-t) 



(5.83) 



The expressions (5.81) and (5.83) for TZ{t) are only con- 
sistent if Ci = and A — 1/3. Hence there are no dust 
solutions in the case when A 7^ 1/3. 



3. A = 0, A : 



Assuming that C2 7^ 0, condition (5.72) reduces to 



nit) + CiTZit) 



Co 



2TV{t) 



which yields 



n[t) 



1/3 



^3C,/. . ^^l/3 ^ ^ 



(^(^-M) 



(5.84) 



(5.85) 



Ci =0. 



The condition that /i be continuous across S reduces 



to 



Thus, 



TZH = c\n 



C2 

7^2 



(5.86) 



H 



7^3 



C2 



2<^+Coe- 
2 
3(to-t)' 



Ci, CiT^O, 
Ci=0, 



(5.87) 



which gives 
a(t): 



aoe-^i*(2CoCi + C2e3C^i*)^ ^i ^ 0, .^ gg. 
ao(io-i)^ Ci=0. ^ ■ ' 



Condition (5.71) yields 






(5.89) 



Condition (5.69) reduces to 



/m(0) 



pr'"itMt))^-PHit) 



n{t) 



CI 



27^G7^(i)5 ■ 



/m(0) 



This equation has no solution for which both pjj and 
Pfj are positive. It follows that there is no physically 
relevant dust solution in this case. 

In summary, in the case A 7^ 1 a static spherical space- 
time cannot be produced by gravitational collapse of a 
homogeneous and isotropic dust fluid. This is different 
from the case A = 1 and that of GR [jj . 



VI. CONCLUSIONS 

In this paper, we have studied gravitational collapse of 
a spherical cloud of fluid with finite radius in the frame- 
work of the nonrelativistic general covariant theory of the 
HL gravity with the project ability condition and an arbi- 
trary coupling constant A. Using distribution theory, we 
have developed the general junction conditions for such a 
collapsing spherical body, with the minimal requirement 
that such junctions should be mathematically meaningful 
in the sense of generalized functions. The general junc- 
tion conditions have been summarized in Table I. Apply- 
ing them to the case where the collapsing star is made 
of a homogeneous and isotropic perfect fluid, while the 
external region is described by a stationary spacetime, 
we have found that the problem reduces to the matching 
of six independent conditions that the Arnowitt-Deser- 
Misner variables {N,N^,gij) and the gauge field A and 
Newtonian prepotential (p must satisfy. For the case of a 
dust fluid, we have shown explicitly that such matching is 
possible only in the case A = Xgr = 1, in which the exter- 
nal spacetimes are described by the Schwarzschild (anti- 
) de Sitter solutions. It is remarkable that the collapse 
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of a homogeneous and isotropic dust to a Schwarzschild 
black hole, studied by Oppenheimer and Snyder in gen- 
eral relativity more than 80 years ago, is a particular 
case. As emphasized previously, our treatments for the 
junction conditions of a collapsing star presented in this 
paper can be easily generalized to other models of the 
Hofava-Lifshitz gravity, or more general to any model of 
high-order derivative gravity theories. 
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Appendix A: Functions (Fs)j, and F^-' , 



The geometric 3-tensors F*-' and F^ . defined 



Eq.(2.16) are given by 
1 



{Fo)ij 
{Fih 
{F2h 



9ij, 



gij R + Ri 



-QijR + 2RRij — 2V(iVj)i? 

+2g,jV^R, 
1 



{F3)ij = —^gijR7nnR"^"^+2RikRj—2W "^ {iRj)k 
-|-V Rij + gijWm'^nR™"', 

{Fijij = --gijR + 3i? i?.y — 3V(,jVj)ii! 
+3g,,V^R^, 

+2RRkiRj — V(iVj) {R"^"Rmn) 
— 2V"V(ii?i?j)„ + gij\7 {R"^"'Rmn) 
+V2 {RR,j) + .9,, V„V„ {RR"'") , 

+ \V^ {R^nR'|) + lg^J^kVl {R^R''^) 

-3VfcV(i {Rj)nR"' ) , 
{Fi^j = \g^AVRf~{V,R){VJR) + 2R,,V^R 



(Fsh 



-2V(iVj)V^i? + 2.g,jV*i?, 



V*i?, 



{y^Rmn) (Vji?™") + 2 (VpR.n) (V'^Rj) 

2v"V(,v2i?^-)„ + 2v„ (i?r„V(,i?j;) 

2V„ {Rm(j^ i)R") — 2V„ lR„i(^iV^^ R"^ 



-.g,jV"V'"V^i?, 



(A.l) 



^S,i) = 1^{(2K + W'^)r^^ -2(2Ki + V^Vk^)R' 
- (2Kg - i?) (2K'^ + V'V^V') } , 



-•^Q- 



■jk 



./2N' 



■ \ ■-/2A^* 



^'P 



^(V,3) 



where 



i {2VfcV(7^'' - VV^' - (V,V,4') .9^^} , 



(A.2) 



rj ^ p{{ 2K'^ + V'V^<^ ]--[2K + V'^] g'- 



(A.3) 



The Fij for the spherical spacetime (3.11 are found to 



be 

{Fo)ij 
(Fih, 
iF2h 



p2i/ 2 



2e 



-2i/ 



2e"4'' 



ee^-" + e^-^ - Sr^iy'' 



+ I2r^ {v'f - 1 



818] 






Ge^-" + e'^" + 4i^'^^'>r^ + 24r^ {i^') 



3^/\3 



(F, 



3)ij 



2rv' {-ie^'' + Ur^iy" + 7) - 7 



4e -f e — or z^ 



^^275 



+ 9r2 (i/')' - 5 






-4i/ 



4e2''-)-e4'' + 3z/(3)r3 + 18r3(zy')' 



(^4). 



-rv' (-46^" + 21r2:/" + lO) - 5 

4e-4'' 
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jj, 



s^'' -1^ 2ru' - 1] 



22e 



2i^ I „4i/ 



24r^i/ 



2,/' 



(5[<5; 



-f40r2 (jy')^ - 2 (e^'' - l) rv' - 23 

+ *^ |240r4 (i.')^ + 4 (1862-^ - 17) r^ {v') 

-12r2 {v'f (-lle^^^ + 22r2i." + 15) 
+irv' \ - Ge^" + 7e'''' + Si/^^V^ 

- 28 (e^^^ - 1) r^v" - l" 

12r4(j.")'-24(e2'^-l)r2i." 
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+ {e^" - 1) (226^'' + e""- + Qu^^\^ - 23) 

+ {e^" - 1) (lee^-^ + e"'' - Wr'^v" - 17) 
+ (21e2'' - 17) r^ {v'f 

~^rv' {-le^" + 9r^iy" + 7) |(5[(S; 



n 



'ij, 



-r^ {i^'f (2262-^ + 101^.(3)^3 _ 102) 

-2r^ {v'f (-Ge^" + IGSr^i/" + 45) - 14 }fl. 



(Fsh 



+ ^^L8r^{^"f + l80r^{i,'f 



r\" 



+ (e^'' - 1) (326^'^ + 2e'*'^ + Iv^'^^'^ - 34) 

+21 (262^^ - 1) r^ {v'f - 28 (e^"^ - l) 
-r^ [v'f (-77e2'' + l%Sr^v" + lOl) 
+rv' [3 (-^e^" + he^" + dv^^'^r^ + 3) 



6 (2e2- + e4- _ ^(4)^4 _ 3^ 

+Abr'^ (y"f + 120r4 {v'f + lOr^ (i/')^ 
-8 (e^- _ 4) r^^y" 

-r^ (i.')' (-12e2'' + 246r2z." + 77) 

+2rv' (86^" + 30i^(3V3 - Sr^i/" - 32) j(S[5J 



-61/ 



(496^'^ - 41) 



^, 



-4iy 



(^6). 



- (e^'' - 1) (1362" + e^'' - Gr^i." - 14) 
-9e2''r2 (i/')' + 6rz/' (-2e2"^ + ^r'^v" + 2) U[<5; 

+ ^|l5r4(i.")' + 150r4(z.')* 

+ (e^'' - 1) (266^'^ + 2e'*'^ + Si/f^V^ - 28) 

+ (iSe^'' + 25) r^ [v'f - 12 {(?'' - l) r 
-3r2 (z.')' (-lie'" + 55r2i/" + 12) 
+3ri/' 



-360r5 (i.')' - 3r4 (v"f - SOr^ {v'f 

_4e2'-,y(3)r3 ^ 16j,(3)r3 

+r^i." flGe^'^ + 75v^^'>r^ - 64) 

+2r3 (i/')3 (-1262" + A89r^v" + 113) 

-r^ {v'f (446^"^ + 303^.(3)^3 _ 33^2^,, _ 269) 

-rv' hsir^ {v"f - 2 {Ue^"" ~ 85) r^v" 
(8e2'' + e4._g^(4)^4_25) 



+6 



36 ^r^ 



+3j,(3)^3 



where v' = dv/dr and fiy = Jf (5^ + sin^ 



(A.4) 



^v" 



Appendix B: Proof of (4.11 )-(4.12 | 



-12e2'-'+4e'"-' + 51/(3)^3 



Let i^ be given by (4.8). We will show that the equa- 



He"' - I) r^v" 



tion i^ = is equivalent to the conditions (4.111 and 



Q. 



'jj, 



-4iy 



(^7), 



(4.12). It is clear that the equation F = is equivalent 



+40r'' {v'f + Arv' he^"" + bv^^K^ - &) 
^2r^v'f {-Se^" + 41r^v" + 15) 
-4 (e^'' - 3) r^z." - 7U[,5; 

-^^ll2e^- + 2e^'' + v^^^r^ 

+12Gr^ {v'f + 262"^ 1/(3)^3 _ g^(3)^3 

+rv' [2462" + e^" ~ Wv'^^^r^ + 127r^ {v"f 

-2 {Te'^" ~ 33) r^v" - 57 

^r'^v" (Se^- + 25^.(3)^3 _ 24) 



to (4.111 together with the condition 



(B.l) 



fe=0 



It remains to show that (B.l) is equivalent to 



Suppose first that (B.l I holds. Then, niultiplyin; 



nt to ( |4.12|. 
Lultiplying( |B.l[ ) 



by $" ^ and using the recursion relation (4.6 ) repeatedly, 
we find 



G'(5(^)($) =0, 0<j< 



(B.2) 



where the function G is defined in a neighborhood of S 
by 

n 

G{x) = ^(-l)'=fc!F^™«(a;)$"-*'-(a;)- (B.3) 

fc=0 



Equation (B.2 1 with j = implies that the restriction of 



G to S vanishes, i.e. Gjs = 0. Equation (B.2) with i — \ 



20 



then gives 
= G<5'($) 



-$<5'($) = --(5($) i.e. - 
$ $ $ 



0. 



In terms of local coordinates {u^} such that u^ — ^ while 
the remaining coordinates {u^}j>2 parametrize the level 
surfaces of $, we have 



= 



G 



dG 



$=0 



Thus, G vanishes to the first order on E. Repeating the 
above procedure n times, we infer that G vanishes to the 
nth order on S: 



G_ 



i.e. 



d^G 
9$J 



0, 



0<j<n. (B.4) 



$=0 



The partial derivatives denoted in the local coor dinate s 
(u^) by jjjj can be expressed invariantly as in (4. 14 1. 



Substituting the expression (B.3) for G into (B.4 1, we 
find 



Vr ^'=^■1 V '^■^^ d^-'^F^-''^^*'^ ^'■^"-fe 



fe=0 



d^i-^ 



a$r 






*=0 

7 \ 9J^("^'=)f-'^™(''^ 



n~kj a$J-(»-fc) 
< j < 71. 



$=0 



Replacing fc by n — fc, we find (4.12 1. 



Conversely, if (4.12 ) holds , the n tracing the above step s 
backwards, we infer that (B.2|, and hence also (B.ll, 
holds. 
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